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Abstract. We propose a new kind of interferometric array that yields images of high dynamic range and large field. The 
numerous individual apertures in this array form a pattern related to a Fresnel zone plate. This array can be used for astrophysical 
imaging over a broad spectral bandwidth spanning from the U.V. (50 nanometers) to the I.R. (20 yum). Due to the long focal 
lengths involved, this instrument requires formation-flying of two space borne vessels. We present the concept and study the 
S/N ratio in different situations, then apply these results to probe the suitability of this concept to detect exoplanets. 
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1. introduction 

Images of high angular resolution and high dynamic range are 
required for the new fields of astrophysics such as exoplanet de- 
tection and cartography of stellar photospheres. Multi-aperture 
interferometry has been used for many years with increasing 
success. However, in the visible domain, interferometric arrays 
are still limited to small numbers of apertures: the maximum 
today being eight, at VLTI. This relatively modest number of 
apertures limits the field-resolution ratio in reconstructed im- 
ages, according to a theoretical limit based on the Shannon the- 
ory of information (Koechlin, Perez 2002). 

To improve the imaging capabilities of interferometric ar- 
rays, we propose a setup allowing the recombination of a very 
large number of beams from very "basic" apertures rectangular 
holes. .The layout of these apertures acts as a diffractive Fresnel 
plate and directly focus the light (combines the beams) into a 
point spread function (PSF) of high dynamic range without the 
need for any reflective or refractive element in the apertures. 
The focal length of such a Fresnel array can vary from 200 
meters to 20 km, depending on the array type and wavelength 
used. This implies space-borne instruments and formation fly- 
ing. 

Imaging diffractive Fresnel plates have been proposed by 
Soret (Soret, 1875) and widely used since (Lipson, Lipson 
& Tannhauser, 1995). Large zone plates have also been pro- 
posed for space borne instruments, using phase or amplitude 
modulation of the incident wavefront to focus light in various 
wavelength domains: submillimetric, I.R., visible (Massonnet, 
2003), X - and gamma-rays (Skinner 2003). 

The interferometric setup proposed here for imaging can 
be regarded either as particular kind of Fresnel plate, or as an 



aperture synthesis array with a very large number of apertures: 
thousands to hundreds of thousands. The apertures are rectan- 
gles punched into a thin metal foil framing the array. The array 
when unfolded in space may have a span of a few to a few hun- 
dred meters (Koechlin 2004). Other interferometer concepts 
involving large numbers of apertures have been developed in 
recent years, such as Carlina (Labeyrie et al, 2004), but with 
classical telescopes forming the individual apertures. 

Fresnel zone plates are chromatic, and so is a Fresnel array. 
For zone plates, the chromaticity issue has been addressed and 
a correction scheme proposed by Falklis & Morris (1989). For 
Fresnel arrays, a similar chromatic correction in the focal in- 
strument provides an achromatic image for spectral bandwidths 
up to ^ = 0.15 as discussed in section Q. The same primary 
array can be used over a wide range of spectral bands whose 
central wavelength is tuned by varying the position of the focal 
satellite along the optical axis of the array. The physical prop- 
erties of the foil defining the primary array limit the observable 
domain to a "global" bandwidth spanning from 50 nm in the 
U.V. to 20 /urn in the I.R.. The limit at short wavelengths is 
the transparency to UV radiation of the thin metal foil defining 
the aperture edges. Towards long wavelengths, the limit is set 
by thermal radiation from the foil itself at the temperature to 
which it can be passively cooled: 40 to 120 K, depending on 
the baffling. 

In the following sections, we present the design of an in- 
terferometric Fresnel array, the percentage of the light going 
into the Point Spread Function, the dynamic range and how it 
can be improved by apodization, the effects of aberrations, the 
achromatizer design and Exoplanet detection simulations. 
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Fig. 1. Example of T c (x,y) for 

ft kmax — 30 Fresnel zones, i.e. 
7200 apertures orthogonal Fresnel array. 

2. Design of the interferometric Fresnel array 

We propose a design where aperture edges follow only two or- 
thogonal directions, the set of apertures (the array) forming a 
large mask. The transmission law T(x, y) of this mask can be 
built as follows. Let us first define functions g and h as: 



8(a) = 1 if 
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and g(a) = otherwise , 
h = \-g 

where a is the distance from optical axis, m is the diffraction or- 
der (m = 1 in our application), k a variable integer: the Fresnel 
zone index and / the desired focal length for the array. 

The function g has an opaque central segment, whereas h 
has a transmissive central segment. If developed radially, g or 
h define circular Fresnel zone plates. In the present case, the 
transmission law T(x, y) of our 2D array is based on an orthog- 
onal development of g and h: 

T c (x,y) = h(x)g(y) + g(x)h(y) 

for a "closed central square" array, and 

T (x,y) = h(x)h(y) + g(x)g(y) 

for the complementary "open central square "array. 

The actual transmission laws depart from T c or T Q , to 
achieve apodization (and mechanical consistency of the grid) 
as described later in the paper. Each aperture (open rectangle) 
in the array can be referenced by its Fresnel zone indices (k x 
and k y ). The term "Fresnel zone" defines an area delimited in 
the aperture plane by two concentric circles. These circles are 



Fig. 2. Computer generated point spread function (PSF) of a 
100-zone (center to edge) 80000-aperture Fresnel array. The 
intensity is displayed at the 1 /2 nd power, to enhance the low 
luminosity regions of the PSF. 

the intersection of the aperture plane with spherical waveplanes 
centered on the focus and whose radii differ by one wavelength. 
The central Fresnel zone is the disc delimited by the smallest 
intersection. The number of zones covered by a Fresnel array 
(as for a filled aperture) corresponds to the number of zones 
crossed from center to edge along a ID line. For a square array 
of size C and k max Fresnel zones, the distances between centers 
of neighboring apertures in the x or y directions are the pseudo 
periods: 
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The apertures cover half the pseudo periods, or less if the array 
is apodized. The focal length of the array is 

/ = C 2 /8k max A 

and the linear PSF half size: 

p — C/Sk max 

For example a 6 m, 600-zone array has an 8 km focal length at 
A = 0.9/j.m, and a linear PSF half size p = 1 .25mm. 

Seen as a diffractive zone plate, this synthesized aperture 
directly forms images of high dynamic range. Most of the light 
that escapes from the central part of the PSF is confined to a 
pair of orthogonal spikes, rather than spread around as for the 
PSF of circular zone-plates. Orthogonal Fresnel arrays are also 
adapted for the "apodized square aperture" approach (Nisenson 
& Papaliolios, 2001), further improving the dynamic range. 

Seen as an interferometric array, this design has the advan- 
tage of not requiring any reflective nor refractive element in the 
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approximated by: <p(x,y) = 2n(y- + y-). 



-PxM 



Fig. 3. evolution of the phase within an individual aperture, as 
seen from common focus. 

primary apertures, and still recombining the beams in a com- 
mon focus. The principle here is to block the part of wavefront 
having non-desired phases as seen from the focus. 

3. Percentage of the incident light going into the 
m th order Point Spread Function 

As Fresnel arrays act by diffraction, only a fraction of inci- 
dent light is sent into the prime focus. As for gratings, there are 
several diffraction orders that get unequal shares, the biggest 
of them being for order zero. We define the "diffraction effi- 
ciency" of a Fresnel array as the proportion of incident light 
going into the first order focus. 

The theoretical diffraction efficiency E clrc for a circular 
Fresnel zone plate when the number of zones tends to infinity 
is: 

Edrc - l/7T 2 m 2 form odd ; E C i rc = for m even and + 

where m is the diffraction order. At the prime focus, which cor- 
responds to the diffraction order m = 1, about 10.1% of the 
light is gathered. 

In the case of the orthogonal array, the diffraction efficiency 
can be derived using Cartesian coordinates: 

The wave contribution at focus can be expressed: 
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(x, y))dxdy. 
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The phase within a given aperture is a quadratic function 
of position, but, when the number of apertures is large, can be 



r p ' /2 C^>' mx my 

¥„, — I I exp(2in — )exp{2in — )dxdy, 

J-nJ2J-l> 



-"Px/2 rPy/2 
-pJ2 J-p Y /2 

which when m + yields 

PxP 



Py 



H 1 



for m odd ; T*,,, = for m even. 



For m odd and + 0, the diffraction efficiency in amplitude for a 
given element of an orthogonal Fresnel array can be expressed 
as: 

vp 4 

and considering a 50% overall "void to total area" transmission 
ratio , the diffraction efficiency in amplitude is 2/jr 2 m 2 . 

Finally, the diffraction efficiency in intensity for an orthog- 
onal Fresnel grid is: 



-'grid 



At order m — 1 , this corresponds to a 4. 1 % efficiency (when 
no apodization is applied). 

The effective luminosity of a diffractive square array of 
size (' is the same as that of a reflective circular mirror of di- 



ameter D — 2C yj -y- = 0.23C on an extended object, and 

D — C V 1 .22 yjEg r jd - 0.497C on a point source. The angular 
resolution and dynamic range of a square array of size C are 
equal or better than that of a circular mirror of diameter D - C. 

4. Dynamic range 

In order to assess the performance of feasible Fresnel arrays i.e. 
with a reasonable number of apertures, we have computed PSFs 
for different test arrays. The dynamic range tests presented in 
this section (see Fig. |4j are a comparison between arrays of 
125, 250, 500 and 1000 zones (respectively 1.2510 5 , 5.10 s , 
2.10 6 and 8.10 6 apertures). Figs.|4]and the following show nor- 
malized diagonal cuts of the PSF They are x-labeled in units of 
resels from the center of the field. A resel, or resolution element 
has an angular extension of A/C. 

We have built these arrays with the transmission law de- 
scribed in section |2 T c (x,y). The PSF of such arrays is the 
square modulus of the Fresnel transform of T c (x,y), noted as 
T c (u, v) in the following. 

As variables x and y are orthogonal, T c (u, v) can be com- 
puted from the functions g and h, which are respectively the 
Fresnel transforms of the g and h one-dimensional transmis- 
sion laws (also defined at section^. 

f c (u, v) = h(u)g(v) + g{u)h(v) 

Relation T c (x, y) = h(x)g(y) + g(x)h(y) still holds for apodized 
arrays, however h + (1 - g) if apodization is applied. 

The fact that two-dimensional arrays can be tested with cal- 
culus reduced to ID for most of the process, greatly improves 
the computation speed and memory requirements, thus allow- 
ing the test of large arrays with little computing power. 
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Fig. 4. Diagonal cut of the PSF, for non-apodized arrays hav- 
ing 125, 250, 500 and 1000 zones: respectively 1.25 105, 5.105, 
2.106 and 8.106 apertures top to bottom curves. The PSF are 
computed by Fresnel transform with a wavefront sampling ad- 
justed to l/25th of the smallest aperture for each array. 



As expected, when the number of zones increases, the nor- 
malized PSF of a square array tends to that of a filled square 
aperture of the same size. The numerical simulations show that, 
in order to achieve a 7.10~ 6 rejection factor (without apodiza- 
tion) at =i eight Airy radii from the center, a k max - 1000- 
zone array is required. A 250 zone-array will only give 7.10~ 5 . 
Increasing k max has drawbacks, such as decreasing the size of 
all including the smallest apertures of the array, at a fixed array 
size. This leads to an increased impact of the element position- 
ing errors on the wavefront quality. Using a complementary 
pattern (T„ instead of T c ) leads to a similar performance in the 
non-apodized case, as one would expect, but causes a modula- 
tion reversal in the outer patterns of the PSF. These differences 
in the outer (low light level) part of the PSF have an impact 
on the dynamic range. They are due to the fact that Babinet's 
Theorem is only an approximation in the case of a limited num- 
ber of Fresnel zones. 

The dynamic range of a given array can be greatly im- 
proved by apodization. We have tested different apodization 
laws (see section^- We have also tested the effects on the PSF 
of wave-plane aberrations caused by aperture mispositioning 
(see section[6}. 

5. Apodization 

A way to improve the dynamic range of a Fresnel array at a 
fixed number of apertures is to apodize. An apodized transmis- 
sion T ap (x,y) is for example: 

T ap (x,y) = T c (x,y)Apod(x,y) 



As the array is orthogonal, the apodization function can be sep- 
arated into x and y: 

Apod(x,y) = Apod x (x)Apod y (y) . 

Apodization can be done in several ways on a Fresnel ar- 
ray. One of them is a transmissive or reflective element at a sec- 
ondary (cooled) pupil plane in the focal setup. Another would 
be to modulate the aperture positions in order to obtain a phase 
effect similar to what is proposed by Guyon (2003) for Phase- 
Induced Amplitude Apodization. 

In this paper we present apodization by modulating the 
apertures areas at the primary array level. Apodizing this way 
somewhat worsens the I.R. noise contribution of the primary 
array by increasing its radiating surface. However, it simplifies 
the design of focal instrumentation and provides a very robust 
means of apodization. Actually, combinations of the different 
apodizations described above may be used and combined with 
coronagraphic devices. 

The complex amplitude contribution of an aperture is not 
proportional to its area, as the waves from a given aperture 
do not all interfere constructively (see Fig. [3}. Thus, to fit an 
apodization law: Apod(x, y), the dimensions of an aperture cen- 
tered at x c and y c must be: 

Px Py 
a{k x ) = — arcsin Apodx(x c ) and b(k y ) = — arcsin Apod y (y c ). 
n n 

To assess the gain in dynamic range and loss in transmission 
for different apodization functions, we have computed the cor- 
responding PSF and global transmission. Pending an algebraic 
derivation of an optimal Apod(x,y) for Fresnel arrays, akin 
to the prolate apodizing functions for filled square apertures, 
(Aime, Soummer and Ferrari, 2002), following Papaliolios and 
Nissenson (2001) and Soummer (2003), we have evaluated 
functions of the form: 

Apodcaiissix) = exp(-jc 2 /jco 2 ) ; xq defined such that transmis- 
sion becomes «o at the limb: 

Xq = -c/2 -y/log flo 
Apod cos (x) - cos[ 2 ^acos(flo)], 

and Apod cos 2(x) = [cos {^acos( -^flo))] , 
where ao is the residual transmission at the edge. Apodizing 
a Fresnel array with Apod cos i(x) leads to damped variations of 
the PSF secondary peaks but a high background level, and low 
overall transmission (see Fig.|5]and Table[Q. 

In order to reach a better rejection factor close to the center 
of the PSF, we have investigated "broader" apodization func- 
tions and apodization starting from T instead of T c arrays. 

Two results compared to previous Apod cos i{x) are 
presented Fig. |6] The apodizations applied on the wave 
contribution are: 

Apod cos i(x) = sin(j(^asin(Apod cos 2(x)))) 

Apod sqr ,(x) = sin(l(-asin{Apod cos i{x))) ) 
Apodtrig(x) = sin(^sin(jApod cos 2(x)J) 

Apod sqn (x) leads to "ringing" in the outer parts of the PSF 
when ao + 0, whereas ao — causes poor transmission. 
Apod lr i g (x) with oq = 0.1, applied to a T„ type grid, leads to 
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Fig. 5. Normalized PSF (diagonal cut) for a 600 zones T c type 
array, non apodized (dashed line) compared to apodized with 
Apod cos i{x) (flo = 0, full line). 
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Fig. 6. Normalized PSF (diagonal cut) with apodization func- 
tions Apod cos 2 (flo = 0), Apodsqn («o = 0), and Apod, r i g (x) 
(ao = 0.1), respectively full, dashed and dotted lines. 



a 2.1% transmission and a rejection factor better than 3 * 10~ 6 
on the diagonal, at 5 resels from center. This last apodization is 
used when generating the PSF for testing exoplanet detection 
in the next sections. 

From Fig. [7J one can see that the dynamic range is high 
in most of the field. In order to achieve that in a whole field, 
a two-step procedure can be used: an exposure at a given ori- 
entation of the PSF, then a second after a 45 degree rotation 
around the optical axis. Each exposure is then split into eight 
sectors: the four 45 degrees centered on the spikes and the four 




Fig. 7. Three-dimensional representation of the PSF resulting 
from Apod, ,i g . Vertical scale is in powers of 10. 

45 degrees centered on the diagonals. A composite image is ob- 
tained, splicing the four best sectors (diagonals) of each image. 

This two-step procedure is not required, except for high dy- 
namic range imaging: the spikes in a single exposure PSF of a 
Fresnel array contain less energy than those caused by a stan- 
dard telescope spider. 

We have measured PSF residuals in a field defined by the 
composite image described above. At less than 4.5 resels from 
the center of the PSF, imaging an exoplanet is not feasible due 
to the poor rejection factor. From 4.5 to 5.5 resels from the 
center, rejection is better than 6 * 10~ 6 in ^ 80% of the field, 
so imaging is partly conceivable. At 5.5 resels or more from 
the center of the PSF, the rejection is everywhere better than 
6 * 10~ 6 . Beyond 8 resels, the rejection gets better than 2 * 10~ 6 . 
Studies are presently being carried out to determine the optimi- 
mum apodization tradeoff for transmission and dynamic range. 
Apodization will also be studied in association with coronog- 
raphy. 

6. Effects of aberrations 

As individual apertures in the array are void rectangles, only 
aperture positioning and dimensioning can affect the wavefront 
before prime focus. 

The effects of aperture mispositioning on the optical path 
difference (OPD), which is a major challenge in other ap- 
proaches of aperture synthesis, are significantly reduced here: a 
mispositioning Ax of an aperture in the plane of the array leads 
to an OPD error 

Ax 

AOPD = A — . 

Px 

The smaller the pseudo period p x , the higher the impact factor 
of mispositioning will be on the OPD. For a given wavefront 
quality, the highest positioning constraints will be put on the 
external elements. 

Large diffractive arrays with long focal lengths have the 
advantage of large pseudo periods, thus providing wavefronts 
which are much more precise than the positioning of the el- 
ementary apertures that gave rise to them. For the long focal 
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Type of Fresnel plate 


circular 


square 


square 


apodization 


no apodization 


apodization 


Value domain of Apod 
from center to edge 




1 -> 


1 -> 0.1 


Applied transmission law 




Apod cos i 


Apod sqn 


Apod lrig 


Grid "positivity" 


opaque central pattern 


void c. p. 


Transmission 


10.1% 


4.1% 


0.25% 


0.90% 


2.1% 



Table 1. Transmission efficiency at focus of diffraction order m — 1 for different types of diffraction patterns and apodization 
laws 





Fig. 8. top: Contours of the PSF resulting from the above ob- 
serving procedure. Shaded areas indicate where the rejection 
factor R is poor: R > 6.10 6 , and white ones indicate where 
R < 6.10~ 6 . The two circles mark 4.5 and 5.5 resels distances 
from the center of the PSF. bottom: Contours and limitations 
of the same PSF combination for a R — 2.1CT 6 threshold (same 
scale as left Fig.). The two circles mark 5.5 and 8 resels from 
center. 



arrays considered here, the OPD errors on a plane wave cross- 
ing a Fresnel array can be approximated as follows (Massonnet 
2003): 

xAx + yAy 



AOPD 



2/ 



Az 



x 2 + y 2 



4/ 2 



where Ax Ay and Az are deviations of an aperture with regard to 
its nominal position, respectively in the plane of the objective 
(x, y) and perpendicular to it (z). 

As the focal length of a Fresnel array is / = C 2 /%nA, where 
C is the side of the square array and n - k max the number of 
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Fig. 9. Diagonal cut of the PSF of an array with "jigsaw" per- 
turbations on the element edges, corresponding to A/4 PTV 
on the wavefront from the smallest (outer) elements. The PSF 
of a perfect array is drawn for comparison. These two curves 
closely overlap. 



Fresnel zones, we have: 



4A 



AOPD = — (k x Ax + k y Ay) + ... 

The coefficients affecting Ax and Az are very small, e.g. 10~ 3 
and 10~ 6 respectively at the edge of an array 6 m in size and of 
focal length 2 km. By comparison the coefficient for Az is 2 in 
the case of a mirror and n - 1 in the case of a refractive material 
of index n. 

We studied the effects of two types of deformations that 
may occur on an array: a random "jigsaw" perturbation of the 
aperture edges and a "parallelogram" shear. 

For "aberrated" arrays, a 2D transmission law is generated 
and arbitrary errors added to the element positioning, then a 
complete 2D Fresnel transform is computed. 

In the jigsaw deformation, the maximum amplitude of per- 
turbations is constant throughout the array: A/4 peak-to-valley 
on the wavefront from the smallest (outer) apertures in the 
array. This corresponds to a cut error of =s 0.3mm on a 6m 
array with 600 Fresnel zones. Due to limits in computing 
equipement, these results were obtained for an array with only 
100 zones (center to edge) i.e. 80000 apertures. 
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Fig. 10. 2-D PSF of an array with a "parallelogram" shear 
distortion. The square root of the PSF intensity is displayed. 
Distortion effects are clearly visible on the diffraction spikes. 
The Strelh ratio is a 71%. 
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Fig. 11. Diagonal cut PSF of an array with a "parallelogram" 
shear distortion (dashed line) compared to the PSF of a per- 
fect array (full line). The central part of the PSF might become 
unusable depending on the target, but the outer peaks are only 
slightly brighter for the disturbed PSF. 



Another possible perturbation in a Fresnel array may be 
caused by a "parallelogram" shear on a non-perfectly rigid 
frame. Numerical simulations were performed on a 100-zone 
array, with a corner shift of one half edge period, resulting to a 
shift of A/2 for the wavefront from the most external patterns. 
The PSF has been computed and we see that this distortion af- 
fects the dynamic range close to the central peak ( Figs.llOland 

EH). 

7. Achromatizer design and local bandwidth 
limitations 

An achromatisation scheme for holographic lenses ( Fig. 1121) 
has been studied and published by Falklis & Morris (1989). 
It uses an achromatic optical element, acting as a field lens, 
which reimages the primary Fresnel zone plate on a secondary 
one in a pupil plane. This secondary zone plate operates at 
order m — -1 and compensates for the chromatic aberration of 
the primary. It is coupled to a converging achromatic element, 
lens or mirror, which makes the beam converge into a final 
focal plane. 

For a Fresnel array, a similar design may be used, where the 
primary array is achromatized using a small circular or orthog- 
onal zone plate. Field optics (L2) in the focal setup reimage the 
array on this secondary zone plate (part of L3), which is carved 
on a reflective surface and blazed for efficiency. The size of the 
secondary is planned to be typically 1 /200 of the primary and 
may operate at order m = -2 rather than -1. The field relay 
optics (L2) consist of a Cassegrain or similar two-mirror com- 
bination. 



Contrary to chromaticism corrections made with dispersive 
materials, this diffractive correction works for all wavelengths 
with no approximation. However, the bandwidth is limited by 
the size of the field optics. 

The chromaticism correction setup also acts as a stop that 
blocks the other orders of diffraction, which add unfocused 
light and destroy the high dynamic range: all orders except 
m = 1 are focused by the field optics and fall at, or close to, the 
center of the secondary Fresnel optics, where a mask blocks the 
light. 

Stray light from objects out of the field, not focused by the 
front array, is focused by the field optics and falls on the edges 
of the pupil plane. As the front array and the stray light are not 
focused exactly in the same plane by the field optic, blocking 
the stray light in all cases requires a narrow opaque margin 
surrounding the front array. This additional band is taken into 
account in the formation flying study by Guidotti (2004). 

For a given position of the prime focal plane, only one 
wavelength (Aq) is strictly in focus at that plane . At wave- 
lengths far from Aq, the half width of the defocused PSF can 
be expressed by 

C(A - Ap) 
P(A) = ^A— ' 

where C is the size of the primary array. 

Field optics of size S at the focal plane correct the chro- 
maticism and yield a diffraction limited PSF, but only for the 
light that can be captured by the field aperture S . For wave- 
lengths too far from Aq, one gets p(A) > S and vignetting effects 
occur, which affect both the transmission and the resolution. 
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With a primary array of size C and a field aperture of size 
S , the bandpass AA centered on A$ is given by: 



AA = 2A 



S - field 
~C 



o- 



where field is the linear non vignetted field at the focus. There 
is a tradeoff between field and bandpass. For example, a 6 m 
array with 50 cm diameter focal optics yields from a maximum 
bandpass at zero field: 

AA _ 1 
~h ~ 6 

to a maximum field a at null bandpass: 

°= 2 

and for a k max = 600 zone array, a field I resolution ratio of: 
Sip = 8k max ^ = 400 

For A outside Ao ± AA/2, the image is still achromatic, but 
the transmission and angular resolution decrease. At the center 
of the field the residual transmission is given by: 



!S A 
trans(A) — trans(Ao) — - 



C A-A Q 



As stated by Falklis and Morris (1989) , F2 and F3 must 
comply with the expression: 



j_ _ j_ 

F 2 D B 



+ 7 fand^ = 7 f + 7 f + ^*[ 7 f + -^---=!-] 

D c F3, D c D D D- r L D A D B F, J 



To explore direct exoplanet detection by standard imaging 
with a Fresnel array, we have simulated several exoplanet situ- 
ations, and for each one we have chosen the smallest array size 
that allows detection with a planet signal /noise > 3 in a maxi- 
mum exposure of five hours. We used the PSF of an undistorted 
apodized array. We present the determined configurations and 
the evolution of the planet signal /noise as a function of A, in 
the spectral band for which the system is resolved angularly. 
The planet signal /noise ratio is derived from the planet/ star 
contrast, the flux received from the star, the zodiacal and exo- 
zodiacal light contribution, the angular resolution of the array 
and the thermal emission from the grid forming the array. 

In the proposed instrumental setup, a baffle protecting the 
Fresnel grid from direct sun prevents specular reflections and 
keeps its temperature down to 60K without active cooling. 
The thermal noise due to emission from the array is computed 
with an emissivity of 0. 1 . This noise is small compared to the 
residual star light at the planet position in all cases presented 
below. The baffle and the related orbital requirements at the 
Sun-Earth Lagrangian points L\ and L2 have been studied by 
Guidotti (2004). 

8. 1 . Star contribution 

The residual power from the star, integrated over a solid angle 
of (2A/C) 2 at the planet's position in the image plane is com- 
puted as follows: 



B„ = 



2hv i 



c 2 [exp^-l] 



In the Fresnel array concept, a satellite supports the Fresnel 
array, another one supports the L2 field lens and focal equip- 
ment. Lens L3 is in fact two lenses: a Fresnel achromatising 
zone plate, (£3,) and an achromatic element (re-imaging part, 
L3). For a perfectly achromatic action of L2 and L3, the use of 
mirrors rather than lenses seems appropriate. Note that F2 is 
linked to L3 dimensions, or vice versa. In order not to have a 
too small 2 nd Fresnel lens, i.e. a manageable size for its edge 
elements, a long enough Dq distance can be set to be virtual by 
the use of a Cassegrain configuration for L2. An optical combi- 
nation allowing sufficient off-axis performances for re-imaging 
the most external patterns is presently under study. 

8. Exoplanet detection 

We propose to use a Fresnel array to detect exoplanets in the 
"imaging" mode, i.e. using a quasi whole field rejection of 
the star light as described in section |4] and a high angular 
resolution image of the planet in that field. Other observing 
modes could be used with this interferometric array, such as 
the "Nuller" approach: any synthetic aperture setup that pro- 
vides achromatic and straight null zones in its PSF is poten- 
tially usable. The Fresnel array can be designed to yield one 
(or two perpendicular) null zones extending from the center of 
the PSF, by shifting the aperture positions by a half period in 
symmetric halves or quadrants of the array. 



P, = e. 



; e B e n(-jl C 2 AvTrsRej 



T e , e e , R e are respectively the photosphere temperature, 
emissivity and stellar linear radius. The star is approximated 
by a uniform disc. 

v is the frequency corresponding to the central wavelength 
used, 

Av frequency range corresponding to the bandpass limited by 
the achromatizer. 

D is the distance of the stellar system, C is the size of the 
Fresnel array, 

7>s,and Rej are respectively the combined transmission 
efficiency of the lens and optical train and the rejection factor 
at the planet's position. 



8.2. Zodiacal and exozodiacal light contribution 

The contribution of both the zodiacal and exozodiacal light are 
integrated over the PSF and over the line of sight. We consider 
the exozodiacal contribution to be the zodiacal light of a solar- 
system at 10 Pc. 

The local zodiacal light is sampled at 90° from the sun, 30° 
from the ecliptic plane. The exozodiacal light is considered at 
the same angles, corresponding to an orbit inclination angle 
i = 60° and a planet at elongation. 
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Fig. 12. Achromatizer setup : Dg is the distance from primary array L\ to field lens L2 placed in the primary focal plane, Dc is 
the distance from L2 to the L3 set, and Do is the distance from L3 set to final focal plane. Distance Da, not shown on this scheme, 
is the distance from the target to the front lens. 



We use the data in "Astrophysical quantities" (Cox & 
Allen, 1999), Table 13.7 for the visible at 0.5 yum, and Fig. 13.1 
and 13.2 for the I.R. at 10.9 /mi. For spectral dependence, we 
modulate these values using a 5800 K blackbody (diffusion in 
the visible) plus a 275 K blackbody (emission in the I.R.) 

The zodiacal and exozodiacal sources are considered as op- 
tically thin and extended objects. The data in Astrophysical 
quantities correspond to zodiacal luminance observed from 
within the ecliptic plane. The exozodiacal luminance integrated 
over the line of sight for an "asterocentric" distance D p — 1 
A.U. is twice as high. Its dependance on D p is approximated 
by I(D P )II{\AU) = D- 2 - 5 . 

Extended light sources contribute to the detection noise by 
their integration over all the PSF, including side lobes. In this 
respect, Fresnel arrays are different from other interferometric 
devices such as imaging or nulling setups using diluted aper- 
tures. 

The PSF of a nulling interferometer tends to zero at some 
locations in the image field, placed at the star's position. It re- 
mains high within broad side lobes designed to cover a large 
part of the planetary system, where a planet is searched for, but 
not imaged. 

An imaging interferometer using diluted apertures has 
broad side lobes too. The noise contribution from extended ob- 
jects is affected by these side lobes, even in the case of pupil 
densification. 

A Fresnel array is a particular case of an imaging interfer- 
ometer whose PSF is close to zero for all the field except a 
small region: the central peak. Its side lobes consist of narrow 
and dim spikes. Both the star and the planet are imaged, but the 
contribution from the star at the planet's position is rejected by 
apodization. For the apodisation law considered here, the en- 
ergy percentage in the central peak of the PSF E pea k, has been 
evaluated by numerical integration to E pea k = 77%, less than 
but comparable to the PSF of a solid aperture. 

"Leakage" in the case of zodiacal light refers to the contri- 
bution from these extended sources integrated over the whole 
PSF, side lobes included. With a diluted aperture array, it cor- 
responds to an integration over a large part of the exoplanetary 
system. With a Fresnel array, as in the case of a single aperture 
instrument, the main contribution of an extended source to a 



given point of the image plane is the integration over a solid 
angle defined by the central peak of the PSF. 

For the different exoplanets, array sizes and wavelengths 
considered here, the angular extension of the PSF peak varies 
from 0.001 to 0.2 planet orbit radius. We have taken into 
account the zodiacal and exozodiacal contributions over the 
whole PSF, but approximated them by extended objects of uni- 
form brightness adjusted to the star-planet distance in each 
case. This leads to the following expression for the total power 
received at focus from zodiacal light in an area corresponding 
to the support of the central peak of the PSF: 

Pzz = 4A 2 (L Z + L EZ )AA Trs 

Lz and Lez are the zodiacal and exozodiacal spectral lumi- 
nances in Wm~ 2 S r -1 fimT 1 

Trs is the transmission efficiency of the lens and optical train. 
8.3. Thermal noise from the array 

The thermal noise contribution from the Fresnel grid seen from 
the focal plane is also computed as a power: 

2hv 3 



c 2 [exp 



hv 
KT. 



1] 



P g = 4A 2 e g B g fillAv 

T g is the grid temperature, e g the grid emissivity; fill the 
filling factor for the grid. 

8.4. Exoplanet contribution 

The power received from the planet at the focal plane is com- 
puted from two contributions, thermal and albedo: 



Pi 



Pth + P n 



All simulated spectra are based on the Planck emission at the 
considered grid, planet and star temperatures. We adapt albe- 
dos and emissivities to each situation, but do not take into ac- 
count the variations of albedo or emissivity within the consid- 
ered spectral bands. We do not take into account absorption in 
the planetary atmosphere. 
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For the thermal emission P,/,, we use the same expression 
as for the star, replacing the temperature, emissivity and radius 
by their planetary counterparts, T p , e p , R p . 

Considering that the planet is at maximum elongation, half 
the disc is illuminated. The reflected part is expressed 

(R e \ 2 7T IR P \ 2 2 
Pal = e e B e AVTI^ — j 2 a P \Dj C trsE Peak 

D p is the star - planet distance, a p and R p are the planet albedo 
and radius, respectively. 

E P eak is the percentage of the energy in the central peak of the 
PSF 



8.5. Signal to noise ratio 

The number of detected photons from each source, star, planet, 
zodiacal, exozodiacal and grid, is computed from the integra- 
tion through bandpass of the power received and the quantum 
efficiency of the detector: 



For all cases we consider 
a 2.1% transmission at focus for the apodized Fresnel array, 
a 50% efficiency for the focal optics train, including the achro- 
matizing zone plate, 

70% quantum efficiency for the detector, 
75% filling factor for the grid. 

The planet signal/noise ratio is computed assuming that 
the noise is the standard deviation of the number of photons 
detected from sources other than the planet. As the number of 
"noise" photons involved is in all cases higher than 10 3 , the 
read noise is considered negligible. 

We have seen previously that for a 600-zone apodized 
Fresnel array, the rejection factor is 

- better than 6 10 -6 in a wide part of the field situated between 
4.5 and 5.5 resels from the center of the PSF, 

- better than 6 10~ 6 in the whole field beyond 5.5 resels, 

- jumps to better than 2 10~ 6 for separations larger than 8 re- 
sels, if the angular resolution leads to such a separation. 

The threshold changes abruptly as we shift from one high point 
of the PSF to the next in rank, depending on the angular po- 
sition in the field. As a result, the signal-to-noise curves pre- 
sented in this paper show a discontinuity at wavelengths corre- 
sponding to a 8 resel planet-star separation. 

We present six situations for which we have computed the 
required 600-zone array size for planet detection at at least 3 
cr around a solar- type star at 10 parsecs, in a maximum in- 
tegration time of 10 hours (two five-hour rotated exposures). 
Note that for wide arrays, the number of zones could easily be 
raised beyond 600, leading to better rejection rates and smaller 
required apertures than the conservative approach taken here. 
These situations are 

- a warm (300K) Jupiter at 1 ua, observed with a 6 m array, 
0.5 m field optics, 10 minutes of integration (Fig.ll3i. 

- a cold (150K) Jupiter at 5 ua, 6 m array, 0.5 m field optics, 
10 hours of integration (Fig.ll4>. 
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Wavelength (micrometers) 

Fig. 13. Warm Jupiter at 1 AU from a solar-type star at 10 
Pc, 6m array, 10 minutes integration, wavelengths from 380nm 
to 650nm. Minimal S/N ratio is better than 4 at the shortest 
wavelengths. 

- a Venus (450K) at 0.7 ua, 15 m array, 1 m field optics, 10 
hours of integration (Fig.ll5>. 

- an Earth (300K) at 1 ua, 40 m array, 3 m field optics, 10 
hours of integration (Fig.ll6>. 

- a cold (150K) Jupiter at 5 ua, same 40 m array, same 3 m 
field optics, 10 hours of integration, allowing detection in the 
I.R. (Fig.E). 

- an Earth (300K) at 1 ua, 120 m array, 3 m field optics, 10 
hours of integration, allowing detection in the I.R. (Fig.ll8>. 

The following curves plot the S/N ratio as a function of 
central wavelength, starting at A = 380nm. The maximal wave- 
length is limited by the constraint that the star-planet separation 
is at least 4.5 resels. A vertical line marks the wavelength cor- 
responding to a 5.5 resels star - planet separation. 

One can see that two spectral domains are favorable: one in 
the visible and one around 10 fim. Although pl "™' ratios are 
less favorable in the visible, imaging exoplanets in the albedo 
dominant region improves detectability in three ways: first by 
enhancement of the rejection factor, due to a wider angular sep- 
aration relative to the resolution, second by a stronger signal 
from the planet, and third by avoiding a dependence on the 
planet's temperature. Notice also that typically for cold Jupiter 
and Earth imaging situations (Fig.[2]& Fig. I16i. the imaging 
wavelength domain (^900nm) is far from the 8 and 5.5 resels 
limitations, so imaging such planetary systems at these wave- 
lengths is possible with smaller arrays, at the cost of increased 
exposure time. 

Limitations due to front grid emission are clearly noticeable 
for the standard Jupiter imaged by a 40 m array, showing differ- 
ent signal-to-noise evolutions beyond 8 fim, depending on the 
grid temperature. This shows that for imaging at wavelengths 
up to 18 /mi, a baffle that keeps the grid temperature =s 60K is 
sufficient to detect a Jupiter-like planet. 

In our simulations, zodiacal and exozodiacal contributions 
to the noise remain below 1%, except in the case of a "cold 
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Fig. 14. Cold Jupiter at 5 AU, 6m array, 10 hours integration, 
wavelengths from 380nm to 3.2/mi. S/N ratio is higher than 4 
at a 900nm. 
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Fig. 15. Venus at 0.7 AU, 15m array, 10 hours integration, 
wavelengths from 380nm to Ll/mi. S/N > 3 throughout the 
observable wavelengths. 
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Fig. 16. Earth at 1 AU, 40m array, 10 hours integration, wave- 
lengths from 380nm to 4.3/zm. S/N > 3 at =s 900nm. 
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Fig. 17. Cold Jupiter at 5 AU, 40m array, 10 hours integration, 
wavelengths from 380nm to 21//m. The S/N in I.R. depends 
upon the grid temperature: 40K to 70K, top to bottom curves, 
but remains > 3 out to 17 /im below 60K. 



Jupiter" between 8 and 12 /im, where the zodiacal light contri- 
bution reaches a maximum of 9% at 10 /im, and in the case of 
an exo-earth at 10 yum, where it reaches 4%. 

As the resolution of an array a few meters wide will not 
allow telluric exoplanet detection at IR wavelengths, front grid 
cooling does not seem to be an issue. Large arrays allow very 
high angular resolution and dynamic range observations even 
deeper in the infrared for specific targets, thus benefit from an 
active cooling to reduce their thermal emission. 

9. Conclusion 

In this paper, we have proposed a concept of a high resolution, 
high dynamic range imager that is able to detect and image exo- 
planets. This can be considered as an interferometer with a very 



large number of apertures. It can only be operated in space and 
requires formation flying. It also requires building precision or- 
ders of magnitude less than a reflective surface or standard in- 
terferometric array of the same angular resolution and dynamic 
range. The light collected over a large area is concentrated on 
smaller classical optics (e.g. l/20th of the array size), to form a 
final image. 

This study on exoplanet detection is just one example of 
the applications of a Fresnel grid, which cover a large frac- 
tion of what can be done in the high angular resolution domain 
with a space telescope of equivalent size, such as imaging com- 
pact objects like rj Car, the envelopes of Be stars or the galactic 
center. As this instrument concept provides large fields (10 6 re- 
sels), a high angular resolution and a high dynamic range, it 
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Fig. 18. Earth at 1 AU, 120m array, 10 hours integration, wave- 
lengths from 380nm to 12.9^m. 

allows imaging of targets from the inner solar system to extra- 
galactic objects. 

Many aspects of this system remain to be optimized before 
a space project can be envisioned, such as the achromatizer and 
a better apodization scheme, for example a PIAAC at the focus 
(Guyon et al, 2005). 

A pre-study of the flying characteristics, pointing perfor- 
mances in terms of speed and precision, and ergol requirements 
for a five-year mission at the L2 Sun-Earth Lagrangian point 
is available in Guidotti (2004). A breadboard project is under 
construction at Observatoire Midi Pyrenees and will serve for 
optical tests in the next two years. 
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